Excitations of amorphous solid helium 
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We present neutron scattering measurements of the dynamic structure factor, S(Q,uj), of amor- 
phous solid helium confined in 47 A pore diameter MCM-41 at pressure 48.6 bar. At low temper- 
ature, T — 0.05 K, we observe S(Q,lj) of the confined quantum amorphous solid plus the bulk 
polycrystalline solid between the MCM-41 powder grains. No liquid-like phonon-roton modes, other 
sharply defined modes at low energy (u < 1.0 meV) or modes unique to a quantum amorphous 
solid that might suggest superfiow are observed. Rather the S(Q,u) of confined amorphous and 
bulk polycrystalline solid appear to be very similar. At higher temperature (T > IK), the amor- 
phous solid in the MCM-41 pores melts to a liquid which has a broad S(Q,cj) peaked near uj ~ 
characteristic of normal liquid 4 He under pressure. Expressions for the S(Q,U)) of amorphous 
and polycrystalline solid helium are presented and compared. In previous measurements of liquid 
4 He confined in MCM-41 at lower pressure the intensity in the liquid roton mode decreases with 
increasing pressure until the roton vanishes at the solidification pressure (38 bars), consistent with 
no roton in the solid observed here. 



I. INTRODUCTION 

The superfluid fraction of liquid 4 He, both bulk and 
4 He in porous media, is traditionally and most accurately 
measured in a torsional oscillator (TO). Below a criti- 
cal temperture T c , the TO frequency increases indicating 
that a fraction of the helium mass has decoupled and 
ceased to rotate with the TO. The effect is denoted a 
nonclassical rotational inertia (NCRI) and the decoupled 
fraction is identified as the superfluid fraction. 

In 2004, Kim and Chan-^ reported a similarly NCRI 
in solid 4 He, in both bulk solid 1 and in solid 4 He confined 
in Vycor. 2 Remarkably, below a T c ~ 200 mK, a small 
fraction of the solid apparently decouples in a TO. The 
NCRI has now been widely reproduced in other labora- 
tories in a variety of sample cells However, the magni- 
tude of the NCRI observed varies between zero and 1.5 %, 
depending on how the sample is prepared, quenched or 
annealed, on the 3 He concentration and on other factors. 
This suggests that the NCRI is associated with defects 
in the solid, dislocations, grain boundaries, amorphous 
regions or other defects^ - — The magnitude and charac- 
ter of the NCRI depends on the oscillator frequency^^ 
and shows effect a 15 i 16 not usually associated with super- 
flow. For example, the NCRI is associated with substan- 
tial elastic energy dissipation in the solid described by 
the Q of the oscillator. Indeed in some cases, the ob- 
served frequency shift, Af, and the disipation, AQ, can 
be quite well described by the real and imaginary parts 
of a common dynamic susceptibility as found in purely 
glassy system o 15 ' 16 . 

Day and Beamis h 17 i 18 and others^ have shown that 
the shear modulus, /i, of solid helium increases at low 
temperature with both a temperature dependence and 
a dependence on 3 He concentration that tracks that of 
the observed NCRI. The increase in /i is attributed to 



the stiffening of the solid as dislocations become pinned 
by 3 He at low temperature. A key question is whether 
the Af can be entirely attributed to elastic behavior or 
whether there is some remaining part that must arise 
from other effects such as superfiow. 

Pursuing this question, Maris 20 and Beamish and 
coworkers 2 ^ have shown that in some TOs there is suf- 
ficient solid helium in the torsional rod to explain the 
observed Af completely in terms of clastic stiffening of 
the solid in the rod. However, in many others it cannoti 2 ^ 
and at least part of the observed Af must have some 
other origin. Similarly, Choi and coworkers 2 ^ 3 - using a 
TO which includes DC rotation have shown that the crit- 
ical velocity depends on the DC velocity as expected for 
genuine superfiow. Thus, while the stiffening of the shear 
modulus below T c and elastic effects can account for the 
Af and AQ in some cases, it cannot account for Af in 
all cases 1 * 2 .^ 22 ^ nor a n effects 

Path Integral Monte Carlo calculations predict that 
the superfluid fraction and Bose-Einstein condensate 
fraction in perfect crystalline solid helium are negligibly 
small. 2 ^2S However, a finite and observable superfluid 
fraction and condensate fraction is predicted 2 ^ in amor- 
phous solid helium. The first one to two layers of helium 
on rough porous media walls are amorphous. Typically, 
the solid in porous media grows from the amorphous lay- 
ers inward toward the center of the pores ^ In Vycor and 
aerogel the tightly bound amorphous layers give way to 
crystalline solid after a few layers so that the solid in the 
interior of the pore is crystalline. However, if the pore 
size is small enough, the solid is amorphous throughout 
the pore, as predicted for classical solids^ Specifically, 
m 47 A pore diameter MCM-41 and 34 A pore diameter 
gelsil, we have shown 3 -! that the entire solid is amorphous 
(no Bragg peaks). Since superfiow in amorphous solid 
4 He is predicted, it is interesting to determine whether 
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the amorphous solid has any low lying modes similar to 
the phonon-roton mode of liquid 4 He, that might sug- 
gest Bose-Einstein condensation (BEC) and superflow, 
or whether it has vibrational modes similar to those of a 
typical polycrystalline solid as observed in classical amor- 
phous solids j22r— This is particularly interesting since 
a low energy mode in solid helium in aerogel has re- 
cently been reported and identified as the origin of local 
superfluidity^ In this context we present measurements 
of the dynamical structure factor (DSF) of amorphous 
solid helium. 



II. EXPERIMENT 

The experiment was performed at the Institut Laue- 
Langevin on the time of flight spectrometer IN5. We used 
an incident neutron wavelength of 5 A, which provided a 
spectrometer energy resolution of 85 p,eV. 

The sample cell was a cylindrical aluminum container 
with an inner diameter of 15 mm and a height of 55 mm. 
The MCM-41 powder sample occupied a volume of 7 cm 3 
corresponding to a height of 40 mm in the cell. The up- 
per part of the cell containing bulk He was masked with 
cadmium. The cell was mounted in a dilution refrigera- 
tor built at the ILL that has a base temperature of 40 
mK. The cell has been filled at T=3.8 K and 80 bar, and 
the solid sample between the grains of the MCM-41 was 
formed by the blocked capillary method. 

In the present 47 A pore diameter MCM-41 powder 
sample, the ratio of the volume in the pores, Vp, to that 
between the grains, Vjg, is^ VpjVic = 0.44. Thus ap- 
proximately 30 % of the helium in the beam is in the 
pores, 70 % between the grains. In a 47 A pore ap- 
proximately 30% of the volume in the pore is occupied 
by the tightly bound amorphous solid layers on the me- 
dia walls. Thus, in the present MCM-41 sample at 1.6 
K when there is liquid in the pores, approximately 80% 
of the helium in the cell is solid (between the grains and 
on the pore walls) and only 20% is liquid (in the pores). 
For this reason, the difference in the scattering intensity 
arising from the solidification of the liquid to an amor- 
phous solid between 1.6 K and 0.05 K is expected to be 
small, as observed in Figs. [5] and [3] below. 

III. RESULTS 

To set the stage, we show the phase diagram of 4 He 
confined in 47 A pore diameter MCM-4 1 31 ' 39 and in 25 
A mean pore diameter gelsi l 37 ' 40 ' 41 in Fig. [TJ At low tem- 
perature and pressure, liquid 4 He in 25 A gelsil (and in 28 
A FSM ) is a superfluid, as shown by Yamamoto et alj* 7 - 
and Taniguchi et al4^, respectively. In the superfluid 
phase, the associated Bose-Einstein condensation (BEC) 
is expected to be connected and continuous across the 
sample providing a continuous phase that enables macro- 
scopic superflow. At higher temperature above the super- 
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FIG. 1: Phase diagram of liquid 4 He in MCM-41 (47 A), gelsil 
(25 A), and FSM (28 A). The superfluid (SF) phase in FSM 
and in gelsil is at the lower left, the blue region. The green 
line is the T c of the SF phase in FSM observed by Taniguchi 
et al.— and red line the T c in gelsil observed by Yamamoto 
et ali^i. The yellow region shows the Localized Bose-Einstein 
condensation (LBEC) region in gelsil. In the LBEC region, 
well defined phonon-roton modes are observe d 38 ' 39 , up to but 
not beyond the pressures and temperatures indicated by the 
blue dashed line in gelsil and MCM-41 which is taken as the 
upper limit of the LBEC region. The red solid line is the 
freezing onset of 4 He in gelsil measured by Shirahama et al.— . 
The black solid and dashed lines show the freezing onset of 
4 He in MCM-41 observed 3 - 9 - from the static structure factor, 
S(Q). The dashed red lines show the phase boundaries in bulk 
4 He. The blue dots at p = 46.8 bar mark the temperatures 
at which the present measurements were made 

fluid phase in porous media the liquid forms a localized 
BEC (LBEC) region in which the BEC is localized to is- 
lands. In the LBEC region the BEC is isolated in patches, 
is not extended, and there is no macroscopic superflow 
across the sampl o 42 ' 43 . The LBEC region lies between 
the superfluid and normal liquid phases as shown in Fig 
□ In 47 A MCM-41 and 25 A gelsil helium solidifies 
at pressures, p > 38 bars. The liquid-solid boundary in 
MCM-41 has been determined by cooling the liquid and 
observing a small reproducible increase in the peak height 
of the static structure factor S(Q) on solidification^ In 
25 A gelsil it has been determined from pressure drop 
on solidification^ - No Bragg peaks were observed in the 
confined solid 31 showing that the solid is amorphous with 
no crystalline regions. In larger pore diameter media, for 
example in 70 A gelsil^ and aerogel^ polycrystalline 
solid 4 He regions are observed. 

Fig.[2]shows the net DSF of helium in the present sam- 
ple cell containing MCM-41 at pressure p = 48.6 bar. 
S(Q,w) at Q = 2.0 A -1 and temperatures T = 0.05 
K and T = 1.6 K is shown. At T = 0.05 K we ex- 
pect amorphous solid helium in the MCM-41 pores and 
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FIG. 2: Net S(Q,ui) at Q = 2.0 A _1 of helium in MCM- 
41 at pressure p = 48.6 bar at T = 0.05 K (amorphous 
solid in the pores) and at T = 1.6 K (liquid in the pores). 
S(Q, uj) includes the scattering from the bulk polycrystalline 
4 He lying between the grains of the MCM-41. At T = 0.05 
K, S(Q,uj) at low uj is approximately proportional to uj 2 
characteristic of scattering from a solid. S(Q, uj) has a peak 
at uj ~ 1.3 meV attributed chiefly to the polycrystalline solid. 
At T — 1.6 K there is liquid in the pores with amorphous 
solid adjacent on the pore walls. S(Q, uj) at 1.6 K shows new 
intensity at low uj characteristic of normal liquid 4 He. 



polycrystalline solid between the grains of the MCM-41 
powder. The inelastic scattering (u> > 0) at T = 0.05 K 
is a sum of S(Q,uj) from the amorphous solid (30 % of 
the helium sample in the beam) and polycrystalline solid 
(70 % of the helium in the beam) . At low uj the intensity 
in S(Q,u)) grows smoothly with uj as expected for scat- 
tering from a polycrystal or an amorphous solid in which 
the density of states g(w) is approximately proportional 
to uj 2 at low uj. There is no indication of any low energy 
roton modes or any layer modes as seen in liquid helium 
in porous media. The roton and layer modes of the liq- 
uid in MCM-41 are observed 3 ^! 4 ^ 4 - 7 - at energies 0.5 - 0.6 
meV and 0.4 - 0.5 meV at 34 bar, respectively. The DSF 
in Fig. [2] peaks at u> ~ 1.3 meV. This is consistent with 
the peak in the density of states (DOS) of phonons in 
bulk polycrystalline helium between the grains observed 
previously . 39 1 46 1 4 7 The large peak at uj = is elastic scat- 
tering, S(Q,uj — 0), from the amorphous solid in the 
pores. The elastic scattering from the polycrystalline 
solid is confined to Bragg peaks and the Q in Fig. [2] has 
been selected to avoid these Bragg peaks. 

In Fig. [2] at T = 1.6 K there is new response at low 
uj in S(Q,uj) not seen at T = 0.05 K. Also at T = 1.6 
K the intensity in S(Q,uj) at higher uj (u> > 1 meV) is 
marginally smaller. At Q = 2.0 A -1 , energy transfers 
uj > 1.5 meV are not observable on IN5 at the present 
incident neutron energy used. From the phase diagram 
Fig. [U we see that the helium in the interior of the MCM- 
41 pores is liquid at p = 48.6 bar and T = 1.6 K. The 



FIG. 3: The net S(Q,u) as in Fig. integrated over Q 
between Q\ = 1.8 A _1 and Q2 = 2.2 A~ 1 as defined in Eq. 
(1). The line labeled "Bose Factor" is the data at T = 0.05 
K multiplied by the Bose factor [ng(w) + 1] for temperature 
T = 1.6 K where ub(uj) is the Bose function. The line displays 
the change in S(Q,uj) between 0.05 K and 1.6 K expected 
from thermal contributions alone. At T = 1.6 K thermal 
contributions to S(Q,uj) are still clearly small. 

helium in the one-two helium layers tightly bound to the 
MCM-41 walls remains amorphous solid at T = 1.6 K and 
higher temperatures. Approximately 30 % of the helium 
in the 47 A MCM-41 pores is in the tightly bound layers 
adjacent to the pore walls which remains solid at T = 1.6 
K. Thus approximately 20% of the total helium sample 
in the beam melts to liquid between T — 0.05 K and 1.6 
K. 

The intensity in S(Q,u) at low uj at T — 1.6 K is at- 
tributed to the normal liquid 4 He in the interior of the 
pores. This S(Q, uj) is similar to that of bulk normal liq- 
uid 4 He at 20 bar a 48 ' 49 and of confined normal 4 He in 
MCM-41 at 34 bars 4 ^ as we discuss further below. That 
is, S(Q,uj) is concentrated near uj= with a tail extend- 
ing out to higher to, as in normal liquidsi 50 ' 51 Particularly, 
there are no sharp modes in S(Q,uj) at either 0.05 K or 
1.6 K characteristic of a Bose-condensed liquid^ 

Fig. [3] shows S(uj) obtained by integrating 
S(Q, uj) over a range of Q values, 

S(oj) = J dQS(Q,uj), (1) 
Qi 

between Q 1 = 1.8 A" 1 and Q 2 = 2.2 A" 1 . The purpose 
of the integration is chiefly to improve the statistical pre- 
cision of the data especially at low uj. As in Fig. [21 the 
additional intensity at low uj at T — 1.6 K is attributed 
to normal liquid 4 He in the interior of the pores. The 
blue line in Fig.[3]shows S(Q, uj) of solid helium observed 
at T=0.05 K (bulk polycrystalline helium between the 
grains and amorphous solid helium in the pores) multi- 
plied by the Bose factor [ns(w) + 1] where Ub(uj) is the 
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Bose function for T = 1.6 K (0.15 meV). The blue line 
represents the change in S(Q,u>) of the solid expected 
when T is increased from 0.05 K to 1.6 K if the energies 
and lifetimes of the modes of the solid remain unchanged. 
Clearly, thermal effects make only a minor contribution 
to the observed change in S(Q,w) between T = 0.05 K 
and T = 1.6 K, even at low to for oj > 0. Thus it is reason- 
able to attribute the observed large change in S(Q, lu) at 
low lu to melting of the amorphous solid in the interior 
of the pores to a liquid. As in Fig. [3J no sharp or well 
defined modes are observed in the integrated S(Q,u>). 



0.2 



I a15 

o 

S 
a> 

Q 



< 

<D 
> 

o 

E 
w 



0.1 



0.05 



P = 48.6 bar 




/ 


T = 0.05 K — ~ 



0.5 



1 

co (meV) 



1.5 



FIG. 4: The net S(uj) as defined in Eq. |T) and shown in 
Fig. (J3]) for amorphous solid in the pores and polycrystalline 
solid between the grains at T = 0.05 K summed over all 
detectors. The sum over all detectors does not show any well- 
defined modes at low E < 1 meV, only an S(oj) expected 
from a density of phonon like modes as given by the sum of 
Eq. (|13p describing a polycrystalline solid between the grains 
and Eq. (120[) an amorphous solid in the pores. 

In Fig. [31 the intensity in the energy range 0.6 < ui < 
1.5 meV is very similar at T = 0.05 K and at T = 1.6 
K. It is marginally smaller at 1.6 K than at 0.05 K but 
has a similar energy dependence. Thus the S(Q,lu) in 
this energy range arising from the amorphous solid (T = 
0.05 K) and the liquid (T = 1.6 K) in the interior of the 
pores appear similar but with marginally larger intensity 
from the amorphous solid. 

Fig.|3]shows the S(lu) data arising from the amophous 
solid in the pores and the polyscrystalline solid between 
the grains at T = 0.05 K summed over all detectors. 
The goal is to show that there are no sharply defined 
modes at ui < 1 meV at accessible Q values. Rather 
S (lj) increases uniformly with lj as expected for scat- 
tering from phonons in a polycrystalline solid. The peak 
in the data summed over all detectors at 1.3 meV is 
broader than the peak in S(Q,uj) at Q —2.0 A _1 shown 
in Fig. 1. The peak at 1.3 meV is expected to arise 
chiefly from the polycrystalline solid helium between the 
grains. Solid helium is a highly anharmonic solid so that 
there will be substantial anharmonic broadening in the 



observed S(Q,u>) and S(w). In the Discussion we com- 
pare S(Q,uj) and S(uj) of anharmonic polycrystalline 
and amorphous solids and provide explicit expressions 
for the S(lu) observed from each. 
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FIG. 5: The net S(Q,u) of liquid 4 He in MCM-41 at T = 
1.2 K and T = 1.6 K interpreted as normal liquid. The lines 
are fits of a damped harmonic oscillator (DHO) function. The 
shape of S(Q,to) is independent of temperature. Only the 
integrated intensity in S(Q,u) increases with T (see inset) 
as more amorphous solid melts to a liquid with increasing 
temperature. S(Q, u) of normal liquid 4 He at higher pressure 
peaks near u — as observed in classical liquid a 50 ' 51 . 
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FIG. 6: The energy and width (FWHM) parameters of the 
DHO function that provides the best fit to the normal liquid 
S(Q,u) shown in Fig. [4] The parameters and therefore 
the shape of the normal liquid S(Q,uj) are independent of 
temperature. Only the magnitude of S(Q, u) increases with 
increasing temperature. 

Fig. [5]shows the difference between S(Q,lu) with nor- 
mal liquid 4 He in the pores (at T = 1.2 K and 1.6 K) and 
amorphous solid in the pores at T = 0.05 K, the latter 
multiplied by the thermal Bose factor for each tempera- 
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ture as discussed in Fig. [5j This difference is interpreted 
as the net scattering from the normal liquid at low u>. 
At higher u>, {u>> 0.8 meV) the S(Q,u>) of the liquid 
and amorphous solid could be quite similar and the net 
S(Q,uj) is not well determined. The net liquid S(Q, uj) is 
a smooth function of ui as expected for normal 4 He. The 
S(Q,ui) also peaks near ui — as observed in classical 
liquid o 50 ' 51 . In contrast, the S(Q,uj) of normal liquid 
4 He at SVP (p « 0) peaks at w~ 0.5 meV. Thus normal 
4 He at higher presssure responds much like a classical 
liquid. The lines in Fig. [5] are fits of a damped harmonic 
oscillator (DHO) function, 



Si(Q,u) 



Zq/tt 



1 — cxp(— hLj/ksT) 
AojujqTq 



(2) 



to the data with energy, ujq, width, 2Tq, and intensity, 
Zq, treated as free fitting parameters. The origin of the 
DHO is discussed in the appendices of RefsJH andl52l and 
the DHO is a standard fitting function in the literature^ 
We found good fits to the data for energies and widths in 
the DHO function that were independent of temperature 
within precision (see Fig. [BJ. Only the weight or inten- 
sity in the DHO increased with increasing temperature 
between 0.8 K and 1.4 K. The weight (strictly the prod- 
uct ujqTqZq which fluctuates less than Zq alone since 
uiq is nearly zero) is shown in the inset of Fig. [5j The 
increase in intensity with temperature was attributed to 
the increase in volume of liquid with temperature as the 
solid melts. The intensity in the liquid S(Q,uj) reaches 
a maximum at T = 1.4 K which is interpreted as the 
temperature when melting of the amorphous solid in the 
pores is complete. At lower pressure^! we also found 
the energy and width of the normal liquid S(Q,u>) was 
independent of temperature over the temperature range 
investigated (1.2 < T < 1.8 K). No physical meaning is 
attributed to the energy or width in terms of modes. The 
DHO is only a convenient representation of S(Q,uj). 

Fig. [7] shows the net elastic scattering from the helium 
in the sample cell at T = 0.05 K and 1.6 K. The top frame 
of Fig. [7] shows S(Q,u) at Q = 2.2 A _1 with the elastic 
peak at to =0. The elastic peak in the top frame arises 
chiefly from the amorphous solid helium in the pores or 
on the grain surfaces. The elastic scattering from the 
polycrystallinc solid between the grains is confined chiefly 
to Bragg peaks which are not seen at Q — 2.2 A -1 . At 
T = 0.05 K the pores are full of amorphous solid. At 
1.6 K there is amorphous helium in the first 1-2 layers 
on the pore walls only, approximately 30 % of the 4 He 
in the pores. The elastic peak is clearly larger at 0.05 K 
indicating more amorphous solid than at 1.6 K. However, 
the peak is only approximately 20 % larger rather than 
a factor of 2-3 as might be expected if only amorphous 
solid contributes to the peak. 

The bottom frame of Fig.[JJshows the elastic scattering 
S(Q, uj) at lu — as a function of Q. In this case Bragg 
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FIG. 7: The net elastic scattering from helium in and be- 
tween the grains of MCM-41 at T = 0.05 K and T = 1.6 
K. At 0.05 K there is amorphous solid throughout the pores, 
at 1.6 K on the pore walls only. Top: S(Q,ui) with elastic 
peak at w =0 attributed to the amorphous solid. Bottom: 
S(Q, uj = 0) vs Q with lines through the amorphous com- 
ponent. Bragg peaks from the bulk hep poly crystalline solid 
between the grains can be seen in S(Q,oj = 0) at Q = 2.02 
A- 1 (1000), 2.15 A" 1 (0002) and 2.29 A" 1 (1011). 



peaks arising from the polycrystalline solid lying between 
the grains are observed at Q = 2.02 A _1 and Q = 2.28 
A^ 1 . The blue and red lines in Fig. [TJare guides to the 
eye through elastic scattering arising from the amorphous 
solid at T = 0.05 K and 1.6 K, respectively. Again, the 
intensity from the amorphous solid is greater at 0.05 K 
when the pores are filled with amorphous solid than at 
1.6 K when there is amorphous solid in 1-2 layers on the 
pore walls only. The difference in intensity is comparable 
to that shown in the top frame of Fig. [7] 

Finally, in Fig. [5] we reproduce from Ref. H3 measure- 
ments of S(Q,uj) at lower pressure where there is liquid 
He in the MCM-41 pores at low temperature, T = 0.4 
K. In Fig. [5] we observe elastic scattering from the amor- 
phous solid layers in the pores at uj = 0, as in the top half 
of Fig.[7Jat p = 48.6 bar. The peak at uj = 0.6 meV is the 
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FIG. 8: S(Q, u) at Q — 2.1 A -1 (the roton Q) versus pressure 
at T = 0.4 K from liquid 4 He in MCM-41 and bulk solid 
helium between the grains (from Ref. [47h . The intensity 
arising from the roton of liquid 4 He in MCM-41 decreases 
with increasing pressure until at 37.8 bar there is little or no 
roton. There is no roton at p — 45 bar. The energy of the 
phonon DOS of the solid increases with increasing pressure. 

roton mode of the Bose condensed liquid at Q = 2.1 A -1 . 
The intensity in the roton mode clearly decreases with 
increasing pressure until there is only a very weak mode 
observed at p = 37.8 bar which is immediately below the 
solidification pressure. No roton mode is observed at p ~ 
45 bar. The decrease in intensity of the phonon-roton 
modes in the liquid with increasing pressure and disap- 
pearance of the mode at the solidification pressure p ~ 
38 bar is consistent with no low-energy, liquid-like mode 
observed here in the amorphous solid at 48.6 bar. The 
broad peak at higher energy u = 1.0 - 1.3 meV in Fig. [5] 
is scattering from the phonons of the polycrystalline solid 
lying between the grains. The energy of the peak arising 
from the phonons increases with increasing pressure as 
expected for phonons. The broad peak at ui ~ 1.2 meV 
at p = 37.8 bar shown in Fig. [8] is similar in shape to the 
peak in S(Q,uj) at uj ~ 1.3 meV observed here at 48.6 
bar, as seen in Figs. Q] and 01 In each case this peak is 
interpreted as arising chiefly from the phonons in the hep 
polycrystalline solid between the grains. 

IV. DISCUSSION 

A. Melting in confinement 

Taniguichi and Susuki^ 3 - have recently reported mea- 
surements of freezing and melting of helium in a 28 
A pore diameter FSM which is quite similar to the 
present MCM-41. They observe melting continuously 
over a wide temperature range, quite different from melt- 



ing of bulk helium, but freezing over a narrower temper- 
ature range. In Fig. 5 above we presented the DSF, 
S(Q,lu), of normal liquid 4 He in the present MCM- 
41. Particularly, on warming, intensity in the liquid 
S(Q, to) is first observed at T= 0.9 K, the intensity grows 
continuously until 1.5 K and saturates to a constant value 
at 1.5 K (see inset of Fig. 5). We interpret this as on- 
set of melting at 0.9 K and melting complete at 1.5 K. 
This increase in intensity of the liquid S(Q,u>) over a 
wide temperature range is consistent with the continu- 
ous melting over a wide temperature range observed by 
Taniguchi and Susuki. They also observe a very small 
volume change on freezing, two orders of magnitude be- 
low the bulk value. This suggests that the molar volumes 
of the amorphous solid and liquid in FSM, and probably 
in the present MCM-41, are nearly identical. 



B. Modes in confined liquid 4 He 

In earlier measurement s 38 ' 39 ' 46 ' 47 we have observed the 
DSF, S(Q,w), of liquid 4 He in the present MCM-41 and 
in 25 A diameter gelsil as a function of pressure and 
temperature. Both phonon-roton and layer modes are 
observed at low temperature. The intensity in the roton 
decreases with increasing pressure, as noted above, and 
with increasing temperature. At p — 34 bar, for exam- 
ple, the intensity in the roton disappears at T = 1.5 K 
and a roton is no longer observed above this tempera- 
ture. The solid circles, squares and triangles in Fig. Q] 
show the maximum temperatures and pressures at which 
P-R modes are observed. At temperatures and pressures 
above the dashed black line through the data points in 
Fig. [TJ well-defined P-R modes are no longer observed. 
Since well defined P-R modes exist when there is Bosc- 
Einstcin condensation, the dashed black line is associated 
with the temperature Tbec at which BEC disappears in 
liquid 4 He confined in these porous media. In 25 A gelsil 
Tbec lies above Tq for superflow as observed in a tor- 
sional oscillator^ 7 - For this reason, the temperature re- 
gion Tc < T < Tbec is identified as a region of localized 
BEC with no superflow across the sample as discussed at 
the beginning of section III. 

Equally interesting, the roton energy decreases with in- 
creasing pressure, decreasing from A = 0.74 meV at SVP 
to approximately A = 0.55 meV at 37.8 ba r 38 ' 46 . The en- 
ergy of a single P-R mode cannot exceed twice the roton 
energ y 54 ' 55 , 2 A. If the single P-R mode energy exceeeds 
2A, it has sufficient energy to spontaneously decay to 
two rotons and the mode will be broadened and not ob- 
servable as well-defined mode. Thus at higher pressure a 
well-defined liquid P-R mode exists at low energy only, 
at wave vectors in the phonon region and in the roton re- 
gion only. For example, at p > 38 bar (at solid pressures) 
a well-defined P-R mode exists in the liquid at energies 
2 A < 1.1 meV only, i.e. at wave vectors in the phonon 
region Q < 0.6 meV and in the roton region 1.7 < Q < 
2.4 A" 1 only. Thus, even if there were trapped liquid in a 
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solid at higher pressure, we would not observe a complete 
liquid P-R mode. 

In summary, our previous finding that the intensity in 
the liquid roton mode decreases with increasing pressure 
and goes to zero at p ~ 38 bar is consistent with the 
absence of a mode at roton energies in the amorphous 
solid in the present measurements. 



C. Amorphous Solid Helium 

In this section we compare expressions for the dynamic 
structure factor (DSF), S(Q,ui), of an anharmonic poly- 
crystalline and an anharmonic amorphous solid. The goal 
is to show that we expect the S(Q,uj) of solid helium in 
these two structures to be very similar when S(Q,lu) is 
integrated over a range of Q values as in Eq. ([1]) . The 
integrated S(uj) may be somewhat more sharply peaked 
as a function of w in a polycrystalline solid because of the 
coupling between the scattering wave vector Q and the 
wave vector, q, of the phonons in the polycrystalline solid. 
This coupling limits the number of phonons that can con- 
tribute to S(ui) in a polycrystal. However, because of (1) 
anharmonic effects which broadening the phonons sub- 
stantially, (2) the large vibrational displacements of the 
atoms which make multimode contributions to S(Q,w) 
large, and (3) because we are considering polycrystals 
rather than single crystals, we expect the difference to 
be small. To illustrate, we write out the single mode 
excitation term, S\(Q,u>), of S(Q,ui) for polycrystalline 
and amorphous solids below. Higher order terms are dis- 
cussed in an Appendix. 

The total S(Q, to) is the sum of an elastic scattering 
term, So(Q), and a series of inelastic terms represent- 
ing excitation of single modes, two modes, interference 
between those modes via anharmonic terms, and higher 
order mode processes, 

S(Q, u) = S Q {Q, 0) + Si(Q, u) + S INT (Q, u) + S 2 (Q, u) - 

(3) 

In a crystal of N atoms where there is periodic trans- 
lational symmetry, we express the vibrational displace- 
ments ui (t) of the atoms I from their lattice points Ri as 
a superposition of waves (phonons) in the crystal of well 
defined wave vector q and polarization index A. With 
this expression, the product [Q.ui(t)] that appears in the 
DSF (see Eqs. (f2"2"j) and (|2l))) in the appendix) is, for one 
atom per unit cell, 



[QMt)} 



^E/^'^^'^AaW- (4) 

qX 



In Eq. dU f q x = (h/2muj q \) 1 / 2 where ui q \ is the fre- 
quency and e q \ is the polarization vector of the wave 
and m is the atomic mass. A q \(t) is the mode annihila- 
tion operator for the phonon (q — 1 to N and A = 1 to 
3). 



In randomly oriented polycrystals of solid helium, we 
assume that, averaged over the polycrystals, 



[Q- e qx] )poly — 



(5) 



We restrict ourselves here to positive ui and low tem- 
perature too » kT so that the number of the thermally 
excited phonons is small. In this case, using the substitu- 
tion (j4]) and the average (0, the one-phonon component 
of 5(Q,w) is 



3 N 



2lTU) q \ 



(6) 



in which I\{Q) is the one-phonon weight factor, I\(Q) = 
d 2 (Q)(HQ 2 /2m), d 2 (Q) is the Debye-Waller factor and 
A(qX, ui) is the one-phonon response function. In an an- 
harmonic crystal such as solid helium, A(qX, ui) has the 



A(q\,Lj) 



8^ 2 A r(gA, W ) 



U 9A 



2w 9A A] 2 



(7) 



where T = T(qX, u>) is the half width of the phonon group 
and A = A(q\,uj) is a shift in the initial frequency ui q \ 
arising from anharmonic effects. In solid helium, the to q \ 
are typically some positive anharmonic frequencies such 
as self-consistent harmonic frequencies. The T(qX,uj) is 
large so that A(qX,uj) is a broad function of w, especially 
for phonons near the Brillouin zone edge o 57 ' 58 . In a har- 
monic limit in which T — > 0, the A(qX,u>) reduces to 



A(X, quj) = 2n [S(u - u q x) - 6(ui - u q \)] 



(8) 



For positive energies, u> > 0, only the first term in 
A(qX,uj) in Eq. © contributes to Si(Q,uj). 

At p=A8.6 bar the one-phonon weight factor I\(Q) = 
d 2 (Q)(hQ 2 /2m) has a broad peak at Q ~ 2.3 A' 1 . Thus 
•for polycrystalline solid helium, we expect Si(Q, ui) to be 
largest at Q values where 1\ (Q) has its maximum and at 
energies where ui q \ is small. 

The phonon density of states (DOS) of a harmonic 
crystal is defined as 



J q x) 



(9) 



qX 



where the sum is over one complete Brillouin zone. The 
St(Q,w) in Eq. ([7|) is proportional to a modified DOS 



3A ^ 2wu} q x 

qX 
qX 



(10) 
(11) 



which is the usual DOS weighted by oj q ^. The gc{u) 
is the modified DOS that is always observed in neutron 
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scattering measurements^. Eq. (| 10[) is the general an- 
harmonic DOS that will be observed in S(Q, uj) of solid 
helium. The second expression, Eq. ([TT]) . for gc(w) holds 
only in a harmonic approximation. Finally, if we inte- 
grate Si(Q,w) over a range of Q values as in Eq. ([I]), we 
obtain 



5i(w) = / dQS x (Q,u) 



(12) 



= [(2irf/Q] I ' dQh(Q)g c (u)S(Q-q-T) (13) 

where we have used NA(Q) = [(2tt) 3 /Q]5(Q) and (fi = 
V/N) is the volume of the unit cell. In Eqs. ([6]) and 
(|13[) , there is a coupling between Q and the phonon wave 
vector q which means that certain values of q only will 
contribute to Sx(uj). A full DOS is not observed. In 
solid helium A(qX,uj) is a broad function in uj, so that 
only an anharmonic DOS given by Eq. (|10j) . substan- 
tially broadened by anharmonic terms, is observed. In 
Fig. 2] we show an S(uS) that is integrated over a lim- 
ited range of Q only. Thus we expect the corresponding 
S'i(w) to be more sharply localized in uj (representing se- 
lected phonons) than the full modified DOS gc{w) given 
by ||TJJ). 

In an amorphous solid, the mean positions Ri do not 
have periodic symmetry. We treat the amorphous solid as 
a large molecule or as a solid that has a single large unit 
cell with all N atoms in the unit cell§£. We express ui (t) 
as a superposition of the normal modes of the molecule 
(numerated by A, A = 1 to 3N) so that Eq. (|4]) becomes, 



IQ-Mt)} 



l 



(14) 



In this case the polarization vectors e\i have an index 
I. We assume that averaged over the N atoms in the 
amorphous solid the polarization vector fulfil, 



IV ' £\l\/Amor — 



(15) 



With this assumption, and restricting ourselves to pos- 
itive uj and low temperatures fku » kT as before, 
Si(Q,u>) for an amorphous solid reduces to, 



S 1 (Q,uj) = S m (Q)I 1 (Q)±^^- 
= S M {Q)h{Q)g A {w) 



A(X,uj) 



(16) 



where Sm(Q) = l/^X); i' exp(— iQ- [Ri — R^]) is a static 
structure factor defined in terms of the mean positions 
Ri of the atoms, h(Q) = d 2 (Q)(hQ 2 /2m) is the one 
mode intensity as before, A(X,uj) is the one mode re- 
sponse function of the form |(7J) and 



9a{u) = 



3N ^ 2ttw a 

A 



A(X,uj) 



3N^lo x 



S(u> - w\). 



(17) 
(18) 



Si(u) = j dQSx{Q,u) 

= §a(u) [ dQS u (Q)h(Q). 



is a modified DOS for the amorphous solid. In a HA, 
A(X,uj) reduces to a S function and Eq. (TT8"]) becomes 
exact. The corresponding Si(u) intergrated over a range 
of Q values is 



(19) 
(20) 



There are two important differences between Si(ui) for 
a polycrystal, Eq. (1131) . and Si(u>) for an amorphous 
solid, Eq. (l20l) . In the amorphous solid we observe the 
modified DOS <m(w) directly in Si(Q,u) unaffected by 
selection of particular modes via a delta function. Also, 
in the amorphous solid Si(Q,u) contains an additional 
factor of Sm(Q), which as h(Q), peaks at Q ~ 2.3 A -1 , 
in solid 4 He at p = 48.6 bar. Thus, we expect Si(Q,u)) 
in the amorphous solid to be somewhat more sharply 
peaked in Q (at Q ~ 2.3 A -1 ) than in the polycrystal. 
If gc(w) an d <7a (w) are similar, we expect the Si(u>) for 
a polycrystal to be somewhat more sharply peaked in 
uj because of the delta function selection of q values in 
the polycrystalline case. However, given that A(X, uj) is 
itself a broad function of uj, we expect the difference to 
be small. 

In addition to the single mode S\(Q,uj) there are 
higher mode terms S 2 (Q,uj), S 3 (Q,uj) in © which are 
not negligible in solid helium. The S*2(Q,w), which is 
discussed in the appendix, is proportional to l2(Q) = 
d 2 (Q)(hQ 2 /2m) 2 and has a broad maximum at some- 
what higher Q values than Si(Q,uj). It peaks in uj at 
higher uj than Si(Q,u>). The chief effect of the higher 
order terms is to further broaden S(Q,uj) and to extend 
S(Q,uj) to higher uj . 

Also, quite generally, at higher uj we expect the 
S(Q,uj) to be independent of phase (liquid or solid) 
and to become quite similar in polycrystalline, amor- 
phous and liquid helium. Ultimately the response at 
high uj arises from high energy interaction between pairs 
of atoms via the hard core of the interatomic partial 
which will be independent of structure and similar in all 
three phases. This high energy region begins at energies 
greater than the collective mode (e.g. phonon) energies. 
Thus we expect S(uj) of amorphous and crystalline 4 He 
to be similar at high uj as well as at lower uj from the 
arguments above. 



V. CONCLUSION 

The dynamical response of amorphous solid helium 
confined in MCM-41 at 48.6 bar as observed in the dy- 
namical structure factor, S(Q, uj), is a smooth function of 
energy (uj) characteristic of a solid that has a vibrational 
density of states approximately proportional to uj 2 at low 
uj. No sharp excitation at low uj similar to the phonon- 
roton mode in Bose-condensed liquid helium is observed. 
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The S(Q, ui) of amorphous and bulk polycrystalline solid 
helium are similar and broad at Q values around 2 A -1 , 
as anticipated for a highly anharmonic solid. Above T 
= 1 K, the amorphous solid melts to normal liquid 4 He. 
The S(Q,u) at Q ~ 2 A -1 of the normal liquid at 
48.6 bar is a broad function that peaks near ui ~ as 
in classical liquids, rather than a broad function peaking 
at a finite ui (e.g. ui = 0.5 meV) as in (more quantum) 
normal liquid He at p ~ 0. 
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VII. APPENDIX 

In this appendix, we present some background on the 
dynamic structure factor, S(Q, u>), to support the expres- 
sions given in section IV B. We begin with the coherent, 
intermediate scattering function, 

S(Q, t) = jr ^2(e- iQ - r ^e iQ ' r i' (0) ) (21) 
1/ 

~ c ;2^J_y^ e -«Q-[Ri-ii i /] e <[Q^iW][Q-« i '(0)]> j 
l,l' 

(22) 

in which r\ (t) = Ri + ui (t) where Ri is the mean posi- 
tion of atom I (lattice vector in a crystal) and ui(t) is 
the displacement of the atom from Ri at time t. Q is the 
wave vector transfer in the scattering and d 2 (Q) is the 
Debye- Waller factor, d 2 {Q) = cxp [— Q 2 (u 2 )/3] in a cu- 
bic crystal. We consider one atom per unit cell and cubic 
symmetry for simplicity. Eq. (1221) omits the anharmonic 
interference terms. The corresponding DSF is 

S(Q,u>) = ^- J dte lult S{Q,t). (23) 

To generate the series Eq. ([3]), we expand the sec- 
ond exponential in Eq. (|22[) in a power series in 
([Q.u l (t)][Q.u l >(0)]) n in the usual way^^. The zero or- 
der term (n = 0) is the elastic term. For a crystal in 
which the Ri has periodic order, the corresponding elas- 
tic DSF is 

S (Q, 0) = d 2 (Q)l e- lQ - [Rl - R i' ] 5{uj) 
l,l' 

= c1 2 {Q)NA(Q-t)6(lo), (24) 

The intensity in Sq(Q, 0) is confined to Bragg peaks at 
the reciprocal lattice vectors r. In the present measure- 
ments involving bulk polycrystalline solid lying between 



the grains of the MCM-41, we generally chose wave vec- 
tors Q to avoid the Bragg peaks. In this way, elastic 
scattering from the polycrystallic solid is not observed. 
Above IK this is not always easy since above IK the poly- 
crystals are continually re-crystallizin g 31 i 61 and a Bragg 
peak may appear during a measurement. 

For an amorphous solid, where the mean positions Ri 
of the atoms do not have periodic order, 

5 (g,o) = rf 2 (g)l^e-^^-« ! ') ( 5M 

1,1' 

= d 2 (Q)S M (Q)6(u J ) (25) 

where Sm(Q) is a static structure factor defined by the 
mean positions, Ri of the atoms. The elastic Sm(Q) is 
not the same as the static structure factor S(Q). We have 
found that the S(Q) of the amorphous solid is somewhat 
more sharply peaked in the peak region than that of the 
liquid but otherwise the S(Q)) of the liquid and amor- 
phous solid are very similar 4i 

The inelastic scattering in which the neutron creates 
or annihilates a single mode arises from the term propor- 
tional to {[Q .ui{t)]{Q .u^ (0)]) (n — 1) in the expansion of 
Eq. (|2"2")l . The corresponding intermediate DSF is, 

Si{Q,t) = d 2 {Q)-j-f ^>2e iQ 'l R '~ R i']{[Q ■ ui (t)} [Q ■ u r (0)] ) . 
l,i' 

(26) 

For a crystal, we expand [Q.ui] in Eq. (|26p in terms of 
phonon modes as in Eq. (|4]). The single mode term of 
S(Q, ui) for a crystal is then, 

Si(Q,co) ~ [n B (u>) + 1] d 2 (Q)l ^/ 9 2 A [Q • e qX ] 2 

x A(q\,ui)NA(Q-q-T). (27) 

where A(q\, ui) is the one phonon response function given 
by Eq. (JT]) in general and by Eq. ((5J) for a harmonic crys- 
tal. Making the assumption Eq. <j5j) for a polycrystal and 
in the limit of low temperatures where ub{u>) may be 
neglected, Eq. (J2TJ) reduces to Eq. ©. 

For an amorphous solid, we expand the [Q.ui] in terms 
of the modes A of the molecule given by Eq. flT4"l) . The 
resulting Si(Q,u) for an amorphous solid is, 

Si(Q,u) ~[n B (to) + l]d 2 (Q)^2e iQ < R '- R ^ 

x i^/l[Q. £A( ][Q'^WA.w). (28) 

A 

With the assumption Eq. (fT5")) for the [Q.e>j] and ne- 
glecting hb(uj), Eq. (|28|) reduces to Eq. (JTTJ) for amor- 
phous solids. As discussed in section IV B, the coherent 
Si (Q, ui) of an amorphous solid is directly proportional to 
the DOS gji (ui) that would be observed in the incoherent 
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DSF. There is no well defined wave vector q and therefore 
no coupling between q and Q as there is in the crystalline 
case. For an anharmonic solid, we expect Si(u) in the 
crystalline and amorphous phases to both be substan- 
tially broadened by anharmonic effects which will make 
the two quite similar. 

The two phonon term S 2 (Q,u>) is the term pro- 
portional to ([Q.u^ (t)] [Q.u^ (0)]) 2 in the expansion of 
Eq. (1221) . Substituting the expansion Eq. §4§ [ n this term, 
we obtain, 

S 2 (Q, w) =<P{Q)±-[n B {u) + 1] E E /"A 2 

qi . A <j2 ,A 

x [Q • ei] 2 [Q • e 2 ] 2 A(Q - q x - q 2 - r)A 2 (12, u) 

(29) 

where 

/CO J ' 
— [n B {J) + l]A{l,J) 
-00 27r 

x [n_e(w — uj ) + 1]A(2,lj — u> ) 

(30) 



and where 1 = qiAi and 2 = q 2 X 2 . For a polycrystalline 
sample in which Eq. (|5|) is assumed to hold and in which 
the temperature is low, S 2 (Q,U)) reduces to 

s 2B .„>= #(0) (g) 2 i(_L) 2 Y. ^ 

x A(12, w)iVA(Q - gi - g 2 - r) (31) 

where 



A(12,u;) = / A(l,J)A(2,u-J) (32) 

J 27T 

Since A(gA, w) is a broad function, the A(12,w) will be 
an even broader function since A(12, uS) is a convolution. 
S 2 (Q,ui) is significant in solid helium and contributes to 
S(Q,ui) at higher uj. When S 2 (Q,uj) is signiheant, we 
expect S(w) given by (fTJ) to be even more similar for 
amorphous and crystalline solids, especially at higher ui. 
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